Stochastic Resonance in a Dipole 
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We show that the dipole, a system usually proposed to model relaxation phenomena, exhibits 
a maximum in the signal-to-noise ratio at a non-zero noise level, thus indicating the appearance of 
stochastic resonance. The phenomenon occurs in two different situations, i.e. when the minimum 
of the potential of the dipole remains fixed in time and when it switches periodically between two 
equilibrium points. We have also found that the signal-to-noise ratio has a maximum for a certain 
value of the amplitude of the oscillating field. 
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The phenomenon of stochastic resonance (SR) may be characterized by the appearance of a maximum in the 

output signal-to-noise ratio (SNR) at an optimized non-zero noise level. Although early studies of SR were restricted 

■ to bistable systems, later developments have shown that this phenomenon is also present in a wider class of situations, 
including threshold devices and fire-and-reset models. In all these cases, a combination of the periodic input signal 
and the presence of an optimized noise may give rise to SR when the system crosses a barrier or a threshold in a 
coherent fashion. This phenomenon has been found in many particular realizations among others a Schmitt trigger 
0, lasers [||, neurons ||l^,|l^ and magnetic particles Q. 

^ In spite of the efforts devoted to show the ubiquity of the phenomenon in many branches of physics, chemistry and 

^ biology, there is an aspect which has not been considered up to now. It refers to the possibility of the appearance 

■ of a maximum in the SNR ]l5| during the relaxation process of a system towards a single minimum, due to the 
f — I combined action of the noise and the periodic signal. In this paper we address this possibility precisely, which opens 

new applications of SR in the domain of the relaxation phenomena. To be more specific, consider a system in an 
equilibrium state, stable or unstable, where no motion is observed. Addition of noise removes the system from this 
situation in such a way that a force field acts on it. If this force varies periodically in time, then SR may occur. As 
a physical realization of the phenomenon we have shown that SR may take place in a single dipole under a periodic 
external field. 

• Let us consider the relaxational dynamics of a single dipole in an external field which may be described by 

the Langevin equation governing the dynamics of the unit vector R along the direction of the dipole moment 



'a: '^^(h{t)xR+m)xR . (i) 



Here the total field h{t) = k{l + asm{u>ot))z, with z a unit vector, consists of a constant plus an oscillating field and 
it is characterized by the parameters k and a. The noise term ^(t) may originate from the presence of an external 
random field or from thermal fluctuations. It is assumed to be Gaussian white noise with zero mean and second 
^ moment {£,i{t)£,j{t + r)) = Ddijd{T), defining the noise level D. 

The potential energy of the dipole is V{9) — —h{t)cos6, where 9 is the angle between its direction and the direction 
of the external field (cos(0) — z ■ R). Two qualitatively different situations may occur, one for small (a < 1) and 
another for high (a > 1) oscillating fields. In the first case, the potential has a minimum which remains fixed in time 
at 9 = (Fig. I^a). In the absence of noise, once the dipole relaxes toward this minimum no motion is observed, 
whereas when noise is added, the motion is modulated by the oscillating field. In the second case (Fig. 0b), the 
position of the minimum switches periodically in time between the equilibrium points 9 = and 9 = tt. In the 
absence of noise, starting from any arbitrary initial condition {9q ^ tt) the dipole relaxes towards 9 = 0, in spite of 
this point becoming stable or unstable periodically. When a small amount of noise acts on the dipole, it may leave 
the unstable equilibrium point toward the stable one. We could easily infer that the output signal has a maximum as 
a function of D but it is not a trivial matter to elucidate the behavior of the SNR, because both signal and noise go 
to zero as D decreases. 

Since we are concerned with the magnitude cos(0), which for a magnetic dipole corresponds to the component of 
the magnetization along the field, the averaged power spectrum is given by 

/•OO /'27t/ujq 

P{io)= I dse''^' I {cos{9t+s)cos{9t))dt . (2) 
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Analitycal studies of the system, using linear response theory (LRT) which is vahd for small fields, show a maximum 
in the signal, i.e. in the susceptibility, as a function of the noise level D. However, the SNR is a monotonically 
decreasing function of D. One finds 

SNR ^L{fi)ka^ , (3) 



obtained from the fluctuation-dissipation theorem |19| and the expression for the complex susceptibility |jl8[ 



X(-) = , (4) 

where T(/i) — fi^L'{^)/2kL{fi) is the relaxation time and L{^) — coth(/i) — is the Langevin function, with 
fi = 2h{t)/D. 

Beyond the domain of applicability of LRT, we have found a range of values of the parameters for which the SNR 
exhibits a maximum for a non-zero noise level. Our result is obtained from numerical simulations by integrating the 
corresponding Langevin equation by means of a standard second-order Runge-Kutta method for stochastic differential 
equations j2^j2^. The SNR, computed numerically through the averaged power spectrum, is defined as 

^»='»"'»-^- 

where AUuo) is the area of the peak above the noise floor and N{luo) is the noise background at the frequency ujq- 
In Fig. we have depicted the SNR for finite amplitudes of the oscillating field. We have observed that for small a 
the SNR is a monotonic decreasing function of the noise, as predicted by LRT. Nevertheless, for higher values of a 
the SNR has a maximum at D ^ (Fig. |^). When the amplitude decreases, this maximum becomes less sharp and 
disappears at lower a in a continuous fashion. 

This behavior could also be obtained analytically in the limit of small frequencies, i.e. when the system relaxes faster 
than the period of the external field. The power spectrum may then be computed from the following approximated 
expression for the correlation function 

{cos{9t)cos{9t,)) « LtLt,{l ~ e^"') + (^1 - e^"' , (6) 

where Lt = L{ii{t)), Xt.r = — jl T{ii{s))^^ds and fit = fJ-it). Additionally, to obtain eq. (||) we have used the previous 
analytical expression of the relaxation time, which differs from the exact value only in a few per cent JlSf . The results 
obtained by using eq. ^ in eq. (H) are depicted in Fig. ||. It is interesting to notice that the SNR defined from 
eq. (^), which has dimensions of the inverse of time, has been written in dimensionless units in such a way that it 
coincides with the definition of eq. (^). The small differences between the analytical and simulation results occur at 
intermediate values of D and come from the approximation in the relaxation time. 

In the previous results we have shown a dependence of the appearance of SR on the amplitude of the input signal. 
This dependence contrasts with the results obtained for the bistable quartic potential, in which the maximum of the 
SNR disappears for sufficiently large amplitudes of the oscillating field. Additionally, numerical simulations show a 
qualitative change of the SNR upon varying the frequency of the signal ujq. The results obtained indicate that when 
increasing luq, for a fixed a, the maximum of the SNR disappears, as illustrated in Fig. ^. 

Although, in Fig. || we have depicted the SNR for some values of a > 1, this behavior is similar to that for a < 1, 
the reason being the fact that when the dipole is able to leave the neighborhood of 6* = the noise is so high that 
it destroys its coherent motion. However, fixed loq for large enough a the amount of noise necessary in order for the 
dipole to leave the unstable equilibrium point does not destroy its coherent motion. This behavior is illustrated in 
Fig. ^, for a = 3, fc = 1 and ujq = 0.1 and is similar to the one found in It shows the time evolution of cos(6') 
for the optimal, lower and greater noise levels. For the previous values of the parameters, the SNR (Fig. ||) exhibits 
a maximum a.t D — 0.033. The addition of noise then increases the SNR up to 15dB over its value at the limit of D 
going to zero. For the optimal value of D the behavior of the dipole looks surprisingly similar to that for the bistable 
potential. 

Recalling our attention to Fig. it also presents another even more interesting feature. For low D the SNR is 
not a monotonic increasing function of the amplitude. This kind of behavior, presented here for the first time, implies 
that in some range of values of the parameters, a small signal is more easily detected than a greater one. The output 
signal, however, is always an increasing function of the amplitude. The maximum of the SNR (Fig. I), as a function 
of a, is due to the fact that the noise at the frequency luq increases faster than the signal when a approaches to 1. 
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For D = 0.01 and luq = 0.01 the SNR has a maximum at a ~ 0.7. In such a case, the SNR may be evahiated by 
approximating the dynamics around the minimum by a two-dimensional Orsntein-Uhlenbeck process ||2C|] . When this 
approximation holds the correlation function reads: 

(cos(..,Ocos(.0)«i^^^%^(^ . (7) 

The SNR has been computed as we did with eq. (||). The corresponding result is represented in Fig. ^. In spite of the 
approximation made to obtain eq. (0), this result is in excellent agreement with the one obtained through simulations 
for the range where this approximation makes sense. It must be pointed out that the SNR, for low amplitudes, follows 
the quadratic law, SNR ex a^. Note that the result given through eq. (0) corresponds to the one for a linear system, 
which shows that nonlinearities do not play any role in the appearance of this phenomenon. What is remarkable is 
the great deal of generality of this result, since it may hold for any system around a minimum. 

In summary, we have shown the existence of SR in a dipole, a system which has been widely analyzed due to its 
importance as a simple relaxation model. It is well-known that the susceptibility of the dipole exhibits a maximum at 
a certain value of the temperature (strength of the noise) [ p4| . In this paper we have also shown that a maximum also 
appears in the SNR. The mechanism responsible for the appearance of SR in this monostable system lies on the role 
played by noise in the relaxational dynamics itself towards the stable equilibrium state. As an important feature, we 
have also found that the SNR depends on the amplitude of the oscillating field and reaches a maximum at a certain 
value of this parameter. This result reveals the existence of an optimal amplitude of the signal in order to be detected. 
Although the dipole is of interest in its own right, our results could also be extended to the situation in which the 
system is around an equilibrium state in a force field whose intensity varies periodically in time. 
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wishes to thank the Generalitat de Catalunya for financial support. 
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FIG. 1. (a) Potential V{9) = —h{t) cos{9) for the maximum of h(t) (continuous line) and for the minimum (dashed line), 
with k = 1 and a = 0.5. (b) Same plot, but with a = 1.5. 

FIG. 2. SNR as a function of the noise level D, obtained through simulations (symbols) and computed from eq. ^ (lines). 
The values of the parameters used here are k = 1, uiq = 0.01, a — 0.5 (squares), a = 0.7 (open circles) and a — 0.9 (filled 
circles) . 



FIG. 3. Contour plot of the SNR as a function of both the noise level D and the amplitude of the oscillating field a, for 
k = 1 and the frequencies: (a) LUo/2n = 0.01, (b) LUo/2n = 0.1 and (c) ujo/2tt = 1. 



FIG. 4. Time evolution of cos(6l) for a = 3, fc = 1, a;o/27r = 0.1 and the noise levels: 73 = 3.3 x 10"^ (top), D = 3.3 x 10"' 
(middle) and D = 3.3 (bottom). 



FIG. 5. SNR for a = 3, A: = 1, loo/2tv = 0.1. 



FIG. 6. SNR as a function of amplitude of the oscillating field a, for k = 1, u>o = 0.01 and D — 0.01. Obtained through 
simulations (symbols) and computed from eq. (Q) (solid line). 
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